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Graphene’s isolation launched explorations of fundamental relativistic physics originating from 
the planar honeycomb lattice arrangement of the carbon atoms, and of potential technological 
applications in nanoscale electronics. Bottom-up fabricated atomically-precise segmented graphene 
nanoribbons, SGNRs, open avenues for studies of electrical transport, coherence, and interference 
effects in metallic, semiconducting, and mixed GNRs, with different edge terminations. Conceptual 
and practical understanding of electric transport through SGNRs is gained through nonequilibrium 
Green’s function (NEGF) conductance calculations and a Dirac continuum model that absorbs the 
valence-to-conductance energy gaps as position-dependent masses, including topological-in-origin 
mass-barriers at the contacts between segments. The continuum model reproduces the NEGF 
results, including optical Dirac Fabry-Perot (FP) equidistant oscillations for massless relativistic 
carriers in metallic armchair SGNRs, and an unequally-spaced FP pattern for mixed armchair-zigzag 
SGNRs where carriers transit from a relativistic (armchair) to a nonrelativistic (zigzag) regime. This 
provides a unifying framework for analysis of coherent transport phenomena and interpretation of 
forthcoming experiments in SGNRs. 


Graphene, a single-atom-thin plane of graphite, has 
been the focus of intensive research endeavors since its 
isolation in 2004P The high degree of interest in this ma¬ 
terial originates from its outstanding electronic, mechan¬ 
ical, and physical properties that result from the planar 
arrangement of the carbon atoms in a honeycomb lat¬ 
tice. Indeed graphene is considered both as a vehicle 
for exploring fundamental relativistic physics, as well as 
a promising material for potential technological applica¬ 
tions in nanoscale electronics and optics P 

However, the absence of an electronic energy gap be¬ 
tween the valence and conduction bands of 2D graphene 
casts doubts on its use in nanoelectronic devices. Nev¬ 
ertheless, theoretical studies had predicted that narrow 
graphene nanoribbons (GNRs) can have a large band 
gap, comparable to silicon (~ 1 eV), depending on the 
ribbon’s width W and edge geometry (as well as pos¬ 
sible doping at controlled positions). Pertinent to our 
work, we note that these predictions were made 3 7 for 
GNRs that have atomically precise armchair edges with 
widths W < 2 nm. Consequently, the most recent ad¬ 
vent and growing availability of bottom-up fabricated 
atomically-precise narrow graphene nanoribbonsj-®" 12 in¬ 
cluding segmented- 13 armchair graphene nanoribbons 
(SaGNRs), opens promising avenues for graphene na¬ 
noelectronics and for detailed explorations of coherent 
electrical transport in nanoribbon-based graphene wires, 
nanoconstrictions, and quantum-point contacts. 

Here we report on the unique apects of transport 
through segmented GNRs obtained from tight-binding 
non-equilibrium Green’s function 1 # (TB-NEGF) calcula¬ 
tions in conjunction with an analysis based on a one¬ 
dimensional (ID) relativistic Dirac model. This model 
is referred to by us as the Dirac-Fabry-Perot (DFP) 
theory (see below for the choice of name). In partic¬ 
ular, it is shown that the valence-to-conduction energy 
gap in armchair GNR (aGNR) segments, as well as the 


barriers at the interfaces between nanoribbon segments, 
can be incorporated in an effective position-dependent 
mass term in the Dirac hamiltonian; the transport so¬ 
lutions associated with this hamiltonian exhibit conduc¬ 
tance patterns comparable to those obtained from the 
microscopic NEGF calculations. For zigzag graphene 
nanoribbon (zGNR) segments, the valence-to-conduction 
energy gap vanishes, and the mass term is consonant 
with excitations corresponding to massive nonrelativis¬ 
tic Schrodinger-type carriers. 

As aforementioned, transport through narrow 
graphene channels — particularly bottom-up fabricated 
and atomically-precise graphene nanoribbon#" 13 — is 
expected to offer ingress to unique behavior of Dirac 
electrons in graphene nanostructures. In particular, the 
wave nature of elementary particles (e.g., electrons and 
photons) is commonly manifested and demonstrated in 
transport processes. Because of an exceptionally high 
electron mobility and a long mean-free path, 1 it has been 
anticipated that graphene devices hold the promise for 
the realization, measurement, and possible utilization of 
fundamental aspects of coherent and ballistic transport 
behavior, which to date have been observed, with 
varying degree s of success, mainly at semiconductor 
interfaces) 15 * 16 ^ quantum point contacts,^ metallic 
wiresf^l and carbon nanotubesP^O Prominent among the 
effects that accompany coherence and ballistic transport 
are conductance quantization (in nanoconstrictions) in 
steps of Go = 2 e 2 /h, which have been found earlier 
for quantum ballistic transport in semiconductor point 
contacts 1 # and metal nanowiresf^However, quantization 
signatures were scarcely observed 20 in GNRs fabricated 
with top-down methods. 

Another manifestation of coherent ballistic transport 
are interference phenomena, reflecting the wave nature 
of the transporting physical object, and associated most 
often with optical (electromagnetic waves, photons) sys- 
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terns or with analogies to such systems (that is, the 
behavior of massless particles, as in graphene sheets). 
Measurements of interference patterns are commonly 
made with the use of interferometers, most familiar 
among them the multi-pass optical Fabry-Perot (OFP) 
interferometer,® 1 where the superposition of all the out¬ 
going light waves, bouncing in a cavity from bounding 
partially reflecting mirrors, yields an oscillatory inten¬ 
sity record (interference pattern) which depends on the 
light wavelength and the distance between the mirrors. 

The quest for demonstration of the particle-wave du¬ 
ality of electrons through measurement of quantum in¬ 
terference phenomena associated with electron transport 
in solid-state devices requires materials and configura¬ 
tions having long mean-free paths. This requirement 
limited early experimental wo rk in the late 1980s to 
semiconductor heterostructuresp® 6 where conductance- 
quantization steps in gate-controlled two-dimensional 
constrictions have been observed. Subsequently inter¬ 
ference patterns for nonrelativistic charge carriers in 
the form of conductance oscillations were observed 22 ^ 
(and interprete d 22 * 23 ! as Fabry-Perot-type phenomena) 
in semiconductor nanowires. The advent of 2D forms 
of carbon allotropes has motivated the study of optical- 
like interference phenomena associated with relativistic 
massless electrons, as in the case of metallic carbon 
nanotubetP®! and graphene 2D p-n junctions.®^ (We note 
that the hallmark of the OFP is that the energy sep¬ 
aration between successive maxima of the interference 
pattern varies as the inverse of the cavity length L.) 

For GNRs with segments of different widths, our inves¬ 
tigations reveal diverse transport modes beyond the OFP 
case, with conductance quantization steps (nGo,n = 
1, 2, 3,...) found only for uniform GNRs. In particular, 
three distinct categories of Fabry-Perot interference pat¬ 
terns are identified: 

1. FP-A: An optical FP pattern corresponding to 
massless graphene electrons exhibiting equal spac¬ 
ing between neighboring peaks. This pattern is 
associated with metallic armchair nanoribbon cen¬ 
tral segments. This category is subdivided fur¬ 
ther to FP-A1 and FP-A2 depending on whether 
a valence-to-conduction gap is absent ( FP-Al , as¬ 
sociated with metallic armchair leads), or present 
( FP-A2 , corresponding to semiconducting armchair 
leads). 

2. FP-B: A massive relativistic FP pattern exhibiting 
a shift in the conduction onset due to the valence- 
to-conduction gap and unequal peak spacings. This 
pattern is associated with semiconducting arm¬ 
chair nanoribbon central segments, irrespective of 
whether the leads are metallic armchair, semicon¬ 
ducting armchair, or zigzag. 

3. FP-C: A massive non-relativistic FP pattern with 
1 /L 2 peak spacings, but with a vanishing valence- 
to-conduction gap, L being the length of the cen¬ 
tral segment. This pattern is the one expected 


for usual semiconductors described by the (nonrel¬ 
ativistic) Schrodinger equation, and it is associated 
with zigzag nanoribbon central segments, irrespec¬ 
tive of whether zigzag or metallic armchair leads 
are used. 

The faithful reproduction of these unique TB-NEGF 
conductance patterns by the DFP theory, including 
mixed armchair-zigzag configurations (where the carriers 
transit from a relativistic to a nonrelativistic regime), 
provides a unifying framework for analysis of coherent 
transport phenomena and for interpretation of experi¬ 
ments targeting fundamental understanding of transport 
in GNRs and the future development of graphene nano¬ 
electronics. 

Results 


Segmented Armchair GNRs: All-metallic. Our 

findings for the case of all-metalli(P3 segmented aGNRs 
(when the number of carbons specifying the width is 
J\f w = 31 + 2, l = 1, 2, 3,...) are presented in Fig. [l] this 
lattice configuration is denoted as AAA (mmm). A uni¬ 
form metallic armchair GNR [see Fig. 0i)L exhibits bal¬ 
listic quantized-conductance steps [see Fig. lja)]. In con¬ 
trast, conductance quantization is absent for a nonuni¬ 
form 3-segment aGNR; see Figs. 0b) -0f). Instead of 
quantized steps, a finite number of oscillations appears, 
whose maxima (at the value of unity) maintain a constant 
energy separation. This behavior is indicative of optical- 
like Fabry-Perot multiple reflections of the DW electron 
within the cavity defined by the two contacts (interfaces 
between the segments of different width). The patterns 
in Figs. 0b) -0f) correspond to the category FP-Al. 

The dependence of the conductance on the width of 
the leads relative to that of the constriction is explored 
by comparing the two junctions (exhibiting sharp lead- 
to-constriction interfaces) depicted in Fig. 011) and Fig. 
Bni). Examination of the TB-NEGF conductances for 
these two segmented aGNRs [Fig. 0c) and Fig. 0e)[ re¬ 
veals that a wider lead [Fig. 0ni)] is associated with 
a stronger confinement (sharp conductance spikes) com¬ 
pared to a narrower lead [Fig.0II)[ (oscillations). In the 
DFP results [Fig.0d) and FigiQf)], which reproduce the 
TB-NEGF results, this trend is accounted for by varying 
the length and height of the mass barriers, which gen¬ 
erates either a weak coupling (closed quantunr-dot-like 
conductance spikes) or a strong coupling (open quantum- 
dot-like Fabry-Perot-type oscillations) to the leads. 

Further insight can be gained by an analysis of the dis¬ 
crete energies associated with the humps of the conduc¬ 
tance oscillations in Fig. 0c) and the resonant spikes in 
Fig. 0 e). Indeed a simplified approximation for the elec¬ 
tron confinement in the continuum model consists in con¬ 
sidering the graphene electrons as being trapped within 
a ID infinite-mass square well (IMSW) of length L i (the 
mass terms are infinite outside the interval L\ and the 
coupling to the leads vanishes). The discrete spectrum 
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FIG. 1. Conductance quantization steps (a) for a uniform metallic armchair nanoribbon (I) contrasted to Fabry- 
Perot oscillations (b-f) for 3-segment all-metallic SaGNRs (II and III). The first two columns (employing red colors) 
display TB-NEGF results. The third column (employing green colors) displays continuum DFP results, which reproduce the 
TB-NEGF results in the middle column. (I-III) Schematics of the nanoribbons employed in the TB-NEGF calculations. A 
3-segment segmented GNR is denoted as Hi'' — Af™ — H\ V . with A/] w (i = 1,2) being the number of carbon atoms specifying 
the width of the ribbon segments. In all configurations, the semi-infinite leads are shown in blue color on the far left and far 
right. Note the steps nGo in (a) reflecting full conductance quantization for the uniform nanoribbon. Instead of steps, the 
SaGNR junctions in (b-f) exhibit Fabry-Perot-type conductance oscillations whose maxima are the places where conductance 
quantization is maintained. The effect of the relative widths of the leads and the constriction on the conductance of the 
all-metallic junctions is shown in (c,e) (TB-NEGF) and in (d,f) (corresponding DFP, respectively), illustrating Fabry-Perot 
oscillations (c,d) for a constriction width close to that of the leads [see schematic (II)] and the development of a sharper 
oscillatory pattern [conductance spikes, see (e) and (f)] for a junction with much wider leads [see (III)]. The DFP approach 
is used to analyze the behavior of the TB-NEGF conductance in the energy range of the IGo step [see (b)[. (IV) Diagram 
of the double mass-barrier used in the DFP method [case of massless Dirac-Weyl electrons with mi = 0 and Mi = 0, Mi 
being the carrier mass in the leads (not shown)]. The double-barrier parameters that reproduce the TB-NEGF results were 
L\ = 58ao, 1/2 = lao, vn 2 v\ = t/3 for (d) and L\ = 60ao, L 2 = 6ao, m 2 Vp = t/ 4 for (f). In the case of the wide leads (III), it is 
worth noting that the DFP theory reproduces the gradual widening of the spikes as a function of increasing energy; naturally 
this trend results from the weakening of the confinement effect due to a stronger coupling to the leads for higher energies. 
The horizontal solid lines at G = IGo (blue online) in (d) and (f) describe the DFP conductance obtained when employing a 
potential V(x ) double barrier [similar in shape to the schematic in (IV)] and the assumption = 0; the result is independent 
of the potential barrier’s heights. In all figures (here and below), when a roman number is placed in the same frame along with 
a letter index, it indicates the corresponding lattice or DFP schematic specified by the roman number, ao = 0.246 nm is the 
graphene lattice constant; t = 2.7 eV is the hopping parameter. 

of the electrons in this case is giverPS by 

E n = \J h 2 Vpkn +M 2 vjr, ( 1 ) 

where the wave numbers k n are solutions of the transcen¬ 
dental equation 

tan(fc n Li) = —hk n /(MvF)- (2) 


In the context of Fig. [TJ M = mi = 0 (massless DW 
electrons), and one finds for the spectrum of the IMSW 
model: 

E n = (n + 1/2)ttTlv f /L 1 , (3) 

with n = 0,1,2,.... Remarkably, the energies associated 
with the TB-NEGF oscillation humps and spikes in Figs. 
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[ljc) and 0c) [or Figs. m and 0f ) in the DFP model] 
follow very closely the above relation. Note the constant 
separation between successive energies, 

E n — E n _i = 2E 0 , n = 1,2,3,..., (4) 

which is twice the energy 

E 0 = nhv F /{2L 1 ) (5) 

of the lowest state. 

As is well known, a constant energy separation of the 
intensity peaks, inversely proportional to the length of 
the resonating cavity [here L\, see Eqs. © and © 
above] is the hallmark of the optical Fabry-Perot, reflect¬ 
ing the linear energy dispersion of the photon in optics 
or a massless DW electron in graphene structures. For 
our purpose, most revealing is the energy offset away 
form zero of the first conductance peak, which equals 
exactly one-half of the constant energy separation be¬ 
tween the peaks. In onedimension, this is the hallmark 
of a massless fermion subject to an infinite- mass-barrier 
confinement pl and it provides ultimate support for our 
introduction of mass barriers at the interfaces of the seg¬ 
mented aGNR. Naturally, in the case of a semiconducting 
segment (see below), this equidistant behavior and 1/2- 
offset of the conductance peaks do not apply; this case is 
accounted for by the Dirac-Fabry-Perot model presented 
in Methods, and it is more general than the optical Fabry- 
Perot theory associated with a photonic cavityPS 

In the nonrelativistic limit, i.e., when hk n <C A4v F , 
one gets 

tan(fc„Li) ~ 0, (6) 

which yields the well known relations k n Li ~ nn and 

E n ~ Mv 2 f + n 2 h 2 n 2 /{2ML 2 ). (7) 

For a massive relativistic electron, as is the case with the 
semiconducting aGNRs in this paper, one has to numeri¬ 
cally solve Eq. ([2]) and then substitute the corresponding 
value of k n in Einstein’s energy relation given by Eq. ©• 
It is worth mentioning here that the inoperativeness 
for one-dimensional cases, due to Klein tunnelingP of 
the electrostatic potential V(x) [see horizontal blue lines 
at G = 1G 0 in Figs. m and0f)] was also noted earlier; 
see the curves for normal incidence (labeled 9 = 0) in 
Fig. 2 of Ref. H3 

Segmented Armchair GNRs: All¬ 

semiconducting. Our results for a 3-segment 
all-semiconducting aGNR are portrayed in Fig. 0 [see 
schematic lattice diagrams in Figs. 01) and §n )] ; this 
lattice configuration is denoted as AAA (sss). A uniform 
semiconducting armchair GNR [see Fig. 01)] exhibits 
ballistic quantized-conductance steps [see Fig. 0a)]. 
In contrast, conductance quantization is absent for a 
nonuniform 3-segment (13 — 7—13) aGNR; see Figs. 0b) 
— 0d). Here, instead of quantized steps, oscillations 
appear as in the case of the all-metallic junctions 


presented earlier in Fig. 0 However, the first oscillation 
appears now at an energy ~ 0.22t, which reflects the 
intrinsic gap A/2 of the semiconducting centra l seg ment 
belonging to the class II of aGNRs, specified 7 '^ by a 
width Af w =31 + 1, l = 1,2,3,.... That the leads are 
semiconducting does not have any major effect. This 
is due to the fact that < Af^, and as a result the 
energy gap m\v F of the central segment is larger than 
the energy gap A iiv F of the semiconducting leads [see 
schematic in Fig.0III)]. 

The armchair GNR case with interchanged widths (i.e., 
7—13 — 7 instead of 13 — 7 — 13) is portrayed in Figs. 
0o) -0f). In this case the energy gap of the semicon¬ 
ducting leads (being the largest) determines the onset of 
the conductance oscillations. It is a testimonial of the 
consistency of our DFP method that it can reproduce 
[see Fig. 0d) and Fig. 0f)] both the 13 — 7—13 and 
7—13 — 7 TB-NEGF conductances; this is achieved with 
very similar sets of parameters taking into consideration 
the central-segment-leads interchange. We note that the 
larger spacing between peaks (and also the smaller num¬ 
ber of peaks) in the 7—13—7 case is due to the smaller 
mass of the central segment (0.166t instead of 0.22t). 

From an inspection of Fig. 0 one can conclude that 
the physics associated with the all-semiconducting AAA 
junction is that of multiple reflections of a massive rela¬ 
tivistic Dirac fermion bouncing back and forth from the 
edges of a particle box created by a double-mass bar¬ 
rier [see the schematic of the double-mass barrier in Fig. 
0111)]. In particular, to a good approximation the ener¬ 
gies of the conductance oscillation peaks are given by the 
IMSW Eq. 0 with Mv 2 F = m x v 2 F = 0.22t (13 - 7 - 13) 
or Mv f = miV F = 0.166t (7 — 13 — 7). In this respect, 
the separation energy between successive peaks in Figs. 
0b), 0c), 0e), and 0f) is not a constant, unlike the 
case of the all-metallic junction. 

The patterns in Figs. 0c) and 0f ) correspond to the 
category FP-B. This generalized oscillations cannot be 
accounted for by the optical Fabry-Perot theory, but they 
are well reproduced by the generalized Dirac-Fabry-Perot 
model introduced by us in the Methods. 

Segmented Armchair GNRs: semiconducting- 
metallic-semiconducting. Our results for a 3-segment 
(7—5—7) semiconducting-metallic-semiconducting aGNR 
are portrayed in Fig. 0 [see schematic lattice diagram in 
Fig. 01)]; this lattice configuration is denoted as AAA 
(sms/ The first FP oscillation in the TB-NEGF conduc¬ 
tance displayed in Fig. 0a) appears at an energy ~ 0.23 1, 
which reflects the intrinsic gap A/2 of the semiconduct¬ 
ing leads (with U\ v = 7). The energy spacing between 
the peaks in Fig. 0a) is constant in agreement with the 
metallic (massless DW electrons) character of the central 
segment with A= 5. The TB-NEGF pattern in Fig. 
0a) corresponds to the Fabry-Perot category FP-A2. As 
seen from Fig. 0b), our generalized Dirac-Fabry-Perot 
theory is again capable of faithfully reproducing this be¬ 
havior. 

A deeper understanding of the AAA (sms) case can be 
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FIG. 2. Conductance quantization steps (a) for a uniform semiconducting armchair nanoribbon (I) contrasted 
to Fabry-Perot oscillations (b-f) of two 3-segment armchair GNRs [(II) and (IV)] with both a semiconducting 
central constriction and semiconducting leads (13-7-13 and 7-13-7). (Ill, V) Schematics of the mass barriers used in 
the DFP modeling, with the dashed line denoting the zero mass. The physics underlying such a junction is that of a massive 
relativistic Dirac fermion impinging upon the junction and performing multiple reflections (above miVp) within a particle 
box defined by the double-mass barrier. (c,e) TB-NEGF conductance as a function of the Fermi energy of the massive Dirac 
electrons in the leads, (d) DFP conductance reproducing [in the energy range of the IGo step, see (b)[ the TB-NEGF result 
in (c). The mass-barrier parameters used in the DFP reproduction were Li = 55ao, miVp = 0.22 1, L 2 = lao, m^Vp = 0.5 1. 
The mass of the electrons in the leads was MiVp = 0.1661. (f) DFP conductance reproducing the TB-NEGF result in (e). 
The parameters used in the DFP reproduction were L\ = 53.6ao, miv% = 0.1661, L 2 = lao, rri 2 Vp = 0.51 1. The mass of the 
electrons in the leads was MiVp = 0.22 1. a 0 = 0.246 nm is the graphene lattice constant; t = 2.7 eV is the hopping parameter. 


gained via an inspection of the density of states (DOS) 
plotted in Fig. |3](c) for the total segmented aGNR (cen¬ 
tral segment plus leads) and in Fig. [3][d) for the the iso¬ 
lated leads. In Fig.[3](c), nine equidistant resonance lines 
are seen. Their energies are close to those resulting from 
the IMSW Eq. Q (with L\ = 60.4ao, see the caption of 
Fig. [3} for a massless DW electron. Out of these nine 
resonances, the first five do not conduct [compare Figs. 
[3](a) and § c)] because their energies are lower than the 
minimum energy (i.e., A/2 = A iiVp ~ 0.231) of the in¬ 
coming electrons in the leads [see the onset of the first 
band (marked by an arrow) in the DOS curve displayed 
inFig.§d)]. 

Segmented Armchair GNRs: Effects of hydro¬ 
gen passivation. As shown in Refs. H] and (5} a detailed 
description of hygrogen passivation requires that the hop¬ 
ping parameters t! for the nearest-neighbor C-C bonds at 
the armchair edges be given by t' = 1.121. Taking this 
modification into account, our results for a 3-segment 
semiconducting-metallic-semiconducting aGNR are por¬ 
trayed in Fig. [4] [see schematic lattice diagram in Fig. 
|4][I)]; this lattice configuration is denoted as ” AAA (sms) 
H-passivation.” The first FP oscillation in the TB-NEGF 
conductance displayed in Fig. |4][a) appears at an energy 
~ 0.281, which reflects the intrinsic gap A/2 of the prop¬ 
erly passivated semiconducting leads (with MY = 7). 


The energy spacing between the peaks in Fig. |4][a) is 
slightly away from being constant in agreement with the 
small mass m\Vp = 0.051 acquired by the central seg¬ 
ment with MY = 5, due to taking t' = 1.12 1. As seen 
from Fig.gb), our generalized Dirac-Fabry-Perot theory 
is again capable of faithfully reproducing this behavior. 

A deeper understanding of the AAA (sms)-H- 
passivation case can be gained via an inspection of the 
DOS plotted in Fig. |4][c) for the total segmented aGNR 
(central segment plus leads) and in Fig. |4](d) for the iso¬ 
lated leads. In Fig. |4][c), eight (almost, but not exacrly, 
equidistant) resonance lines are seen. Their energies are 
close to those resulting from the IMSW Eq. (© (with 
L\ = 59.5ao and m\v 2 F = 0.051; see the caption of Fig. 
[ 4 ]) for a Dirac electron with a small mass. Out of these 
eight resonances, the first six do not conduct [compare 
Figs.[4](a) and[4|c)] because their energies are lower than 
the minimum energy (i.e., A/2 = Mliv F ~ 0.281) of the 
incoming electrons in the leads [see the onset of the first 
band (marked by an arrow) in the DOS curve displayed 
in Fig. |4][d)]. From the above we conclude that hydro¬ 
gen passivation of the aGNR resulted in a small shift of 
the location of the states, and opening of a small gap for 
the central metallic narrower (with a width of NY = 5) 
segment, but did not modify the conductance record in 
any qualitative way. Moreover, the passivation effect can 
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FIG. 3. Conductance for a 3-segment nanoribbon with 
a metallic (A/” 2 * = 5) central constriction and semicon¬ 
ducting leads (A/i* = 7). The semi-infinite leads (in blue 
on the far left and far right) are an extension of the left and 
right semiconducting segments with A') 1 ’ 1 ’ = 7; see schematic 
lattice diagram in (I). (II) Schematics of the mass barriers 
used in the DFP modeling, with the dashed line denoting the 
zero mass. The physics underlying such a junction is that of 
a massive relativistic Dirac fermion impinging upon the junc¬ 
tion, which loses its mass upon tunneling in the central seg¬ 
ment and performs multiple reflections within a particle box 
defined by the double-mass barrier, (a) TB-NEGF conduc¬ 
tance as a function of the Fermi energy of the massive Dirac 
electrons in the leads, (b) DFP conductance reproducing the 
TB-NEGF result in (a). The mass-barrier parameters used in 
the DFP reproduction were L\ = 60.4ao, mi = 0, L 2 = lao, 
rri 2 Vp = 0.37 1. The mass of the electrons in the leads was 
A4iVp = 0.23 1. (c)-(d) The total DOS of the junction and 
the density of states in the isolated leads, respectively, ac¬ 
cording to the TB-NEGF calculations. The arrows indicate 
the onset of the electronic bands in the leads. Note that the 
DOS in (c) reveal the existence of five sharp electronic states 
below the onset (at 0.23f = MiVp ) of the first band in the 
leads [see (d)], which consequently do not generate any con¬ 
ductance resonances [see (a) and (b)]. Note further in (c) 
the equal energy spacing between the vertical lines [the five 
solid (red) and four dashed (black) ones] associated with the 
resonances of a massless electron confined within the central 
metallic aGNR segment, a 0 = 0.246 nm is the graphene lat¬ 
tice constant; t = 2.7 eV is the hopping parameter. 


FIG. 4. H-passivation effects in the conductance 
of a 3-segment armchair nanoribbon with a metallic 
(A/ 2 * = 5) central constriction and semiconducting 
leads (Ai W = 7); see schematic lattice diagram in (I). Note 
that the nearest-neighbor C-C bonds at the armchair edges 
(thick red and blue lines) have hopping parameters t! = l.l2t. 
(II) Schematics of the position-dependent mass field used in 
the DFP modeling, with the dashed line denoting the zero 
mass. The physics underlying such a junction is that of a 
massive relativistic Dirac fermion impinging upon the junc¬ 
tion, which reduces its mass close to zero upon tunneling in 
the central segment and performs multiple reflections within 
a particle box defined by the double-mass barrier, (a) TB- 
NEGF conductance as a function of the Fermi energy, (b) 
DFP conductance reproducing the TB-NEGF result in (a). 
The mass parameters used in the DFP reproduction were 
Li = 59.5ao, miVp = 0.05 1, L 2 = 1.5ao, m, 2 V 2 p = 0.30T The 
mass of the electrons in the leads was A4 iv% = 0.28 1. (c)-(d) 
The total DOS of the junction and the density of states in the 
isolated leads, respectively, according to the TB-NEGF calcu¬ 
lations. The arrows indicate the onset of the electronic bands 
in the leads; note the shifts from 0.23f to 0.28t and from 0.42t 
to 0.38f for the onsets of the first and second bands, respec¬ 
tively, compared to the case with t' = t in Fig. [3][d). Com¬ 
pared to Fig. [3] the subtle modifications of mass parameters 
brought about by having t' = 1.12 1 result in having six sharp 
electronic states [see (c)[ below the onset (at 0.28f = Miv%) 
of the first band in the leads [see (d)[, which consequently do 
not generate any conductance resonances [see (a) and (b)[. In 
addition, within the energy range (0 to 0.4t) shown in (a) and 
(b), there are now only two conducting resonances, instead of 
three compared to Figs. §a) and |3][b). a 0 = 0.246 nm is the 
graphene lattice constant; t = 2.7 eV is the graphene hopping 
parameter. 
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FIG. 5. Conductance for ZZZ (all-zigzag, left column) 
and AZA (armchair-zigzag-armchair, right column) 
segmented nanoribbon junctions. See corresponding lat¬ 
tice diagrams in (I) and (II). The 3-segment GNRs are de¬ 
noted as A/r ~ A/ 2 W — A/"i W , with N™ (i = 1,2) being the 
number of carbon atoms specifying the width of the ribbon 
segments. The armchair leads in the AZA junction are metal¬ 
lic (A/j w = 23, class III aGNR). (a)-(b) TB-NEGF conduc¬ 
tance for the ZZZ and AZA junction, respectively, (c) DFP 
conductance reproducing the TB-NEGF result in (a) for the 
ZZZ junction, (d) DFP conductance reproducing the TB- 
NEGF result in (b) for the AZA junction. In spite of the 
different edge morphology, the Fabry-Perot patterns in (a) 
and (b) are very similar. The central zigzag segment con¬ 
trols the Fabry-Perot patterns. According to the continuum 
DFP analysis, the physics underlying such patterns is that of a 
massive nonrelativistic Schrodinger fermionic carrier perform¬ 
ing multiple reflections within a cavity defined by a double¬ 
mass barrier [see diagram in (III)], but with the additional 
feature that V\ = —m\V% and V) = —Miv% are also con¬ 
sidered for segments or leads with zigzag edge terminations 
(see text for details). The mass and Vi parameters used in 
the DFP calculations were Li = 30ao, miVp = 2.23 1 — cEt, 
with c = 7.3, Vl = —miVp, L 2 = l.lao, m 2 Vp = 0.38 t, 
V 2 = —m 2 Vp/ 3, MiVp = 2.30£, Vi = —Miv% in (c) and 
L\ = 29.1ao, miVp = 2.65 t — cEt, with c = 8.4, V% = —miVp, 
L 2 = l.Oao, m 2 Vp = 0.30t, V 2 = — m 2 Vp , Mi = 0, Vi = 0 
in (d). Mi and Vi denote parameters of the leads. E is the 
energy in units of t. ao = 0.246 nm is the graphene lattice 
constant; t = 2.7 eV is the hopping parameter. 


be faithfully captured by the Dirac FP model by a small 
readjustment of the model parameters. 

All-zigzag segmented GNRs. It is interesting to in¬ 
vestigate the sensitivity of the interference features on the 
edge morphology. We show in this section that the rel¬ 
ativistic transport treatment applied to segmented arm- 
chaie GNRs does not maintain for the case of a nanorib¬ 


bon segment with zigzag edge terminations. In fact 
zigzag GNR (zGNR) segments exhibit properties akin 
to the well-known transport in usual semiconductors, 
i.e., their excitations are governed by the nonrelativis¬ 
tic Schrodinger equation. 

Before discussing segmented GNRs with zigzag edge 
terminations, we remark that such GNRs with uniform 
width exhibit stepwise quantization of the conductance, 
similar to the case of a uniform metallic armchair-edge- 
terminated GNR [see Fig. |2ja)]. 

In Fig. §a), we display the conductance in a three- 
segment junction [see lattice schematic in Fig. [5^1)] when 
all three segments have zigzag edge terminations (de¬ 
noted as ZZZ). The main finding is that the central seg¬ 
ment behaves again as a resonant cavity that yields an 
oscillatory conductance pattern where the peak spacings 
are unequal [Fig. [5ja)]. This feature, which deviates 
from the optical Fabry-Perot behavior, appeared also in 
the DFP patterns for a three-segment armchair junction 
whose central segment was semiconducting, albeit with a 
different dependence on L [see Figs. [ 2 ] and Eq. §]. More¬ 
over, from a set of systematic calculations (not shown) 
employing different lengths and widths, we found that the 
energy of the resonant levels in zGNR segments varies on 
the average as ~ ( n/L ) 2 , where the integer n counts the 
resonances and L indicates the length of the central seg¬ 
ment. However, a determining difference with the arm¬ 
chair GNR case in Fig. [2] is the vanishing of the valence- 
to-conductance gap in the zigzag case of Fig. [5ja). It is 
well known that the above features are associated with 
resonant transport of electronic excitations that obey the 
second-order nonrelativistic Schrodinger equation. 

Naturally, one could formulate a continuum transport 
theory based on transfer matrices (see Methods) that use 
the ID Schrodinger equation instead of the generalized 
Dirac Eq. ([8]). Such a Schrodinger-equation continuum 
approach, however, is unable to describe mixed armchair- 
zigzag interfaces (see below), where the electron tran¬ 
sits between two extreme regimes, i.e., an ultrarelativis- 
tic (i.e., including the limit of vanishing carrier mass) 
Dirac regime (armchair segment) and a nonrelativistic 
Schrodinger regime (zigzag segment). We have thus been 
led to adopt the same Dirac-type transfer-matrix ap¬ 
proach as with the armchair GNRs, but with nonvan¬ 
ishing potentials V = r/Mvp, with rj = Q(—E) — O(E), 
where O(E) is the Heaviside step function. This amounts 
to shifts (in opposite senses) of the energy scales for par¬ 
ticle and hole excitations, respectively, and it yields the 
desired vanishing value for the valence-to-conduction gap 
of a zigzag GNR. 

The calculated DFP conductance that reproduces well 
the TB-NEGF result for the ZZZ junction [Fig. [5ja)] is 
displayed in Fig. [5j[c); the parameters used in the DFP 
calculation are given in the caption of Fig. [5j We note 
that the carrier mass (mi) in the central zigzag segment 
exhibits an energy dependence. This is similar to a well 
known effect (due to nonparabolicity in the E — k disper¬ 
sion) in the transport theory of usual semiconductors 












































We further note that the average mass associated with a 
zigzag segment is an order of magnitude larger than that 
found for semiconducting armchair segments of similar 
width (see captions in Fig.[2j and this yields energy levels 
~ (n/L) 2 close to the nonrelativistic limit [see Eq. §]. 
We note that FP pattern of the ZZZ junction belongs to 
the category FP-C. 

Mixed armchair-zigzag-armchair segmented 
GNRs. Fig. [5] (right column) presents an example of a 
mixed armchair-zigzag-armchair (AZA) junction, where 
the central segment has again zigzag edge terminations 
[see lattice schematic in Fig. |5](II)] . The corresponding 
TB-NEGF conductance is displayed in Fig. [5](b) . In spite 
of the different morphology of the edges between the leads 
(armchair) and the central segment (zigzag), the conduc¬ 
tance profile of the AZA junction [Fig. [5|(b)] is very sim¬ 
ilar to that of the ZZZ junction [Fig. 0a)]. This means 
that the characteristics of the transport are determined 
mainly by the central segment, with the left and right 
leads, whether zigzag or armchair, acting as reservoirs 
supplying the impinging electrons. 

The DFP result reproducing the TB-NEGF conduc¬ 
tance is displayed in Fig. 0d), and the parameters used 
are given in the caption. We stress that the mixed AZA 
junction represents a rather unusual physical regime, 
where an ultrarelativistic Dirac-Weyl massless electron 
(due to the metallic armchair GNRs in the leads) tran¬ 
sits to a nonrelativistic massive Schrodinger one in the 
central segment. We note that FP pattern of the AZA 
junction belongs to the category FP-C. 

Mixed zigzag-armchair-zigzag segmented 
GNRs. Finally Fig. [6] presents an example of a mixed 
zigzag-armchair-zigzag (ZAZ) junction, where the cen¬ 
tral segment corresponds to a metallic armchair GNR 
[see lattice schematic in Fig. 01)]. The corresponding 
TB-NEGF conductance is displayed in Fig. 0a)- In 
spite of the different morphology of the edges between 
the leads (zigzag) and the central segment (armchair), 
the conductance profile of the ZAZ junction [Fig. 0a)] 
is controlled by determined the central segment, with 
the left and right leads acting as reservoirs supplying the 
impinging electrons. Thus the conductance peaks are 
close to being equidistant, and the FP pattern belongs 
to the categoty FP-A1. The DFP result reproducing 
the TB-NEGF conductance is displayed in Fig. 0b), 
and the parameters used are given in the caption. 

In a reverse sense compared to the AZA junction 
above, the mixed ZAZ junction here represents also a 
rather unique physical regime, where a nonrelativistic 
massive Schrodinger electron (due to the zigzag GNRs 
in the leads) transits to an ultrarelativistic Dirac-Weyl 
massless one in the central segment. The decisive ad¬ 
vances brought forward by our DFP ID continuum the¬ 
ory can be clearly appreciated by its ability to describe 
the corresponding TB-NEGF results for these highly non¬ 
trivial ZAZ and AZA junctions [compare Figs. 0a) and 
0b), as well as Figs. 0b) and0d)]. Contimuum 2D for¬ 
mulations are unable to describe the important all-zigzag 
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FIG. 6. Conductance for a ZAZ (zigzag-armchair- 
zigzag) segmented nanoribbon junction. See corre¬ 
sponding lattice diagram in (I). The 3-segment GNRs are de¬ 
noted as A/”i M — A4 W — A , with N™ (i = 1,2) being the 
number of carbon atoms specifying the width of the ribbon 
segments. The central armchair segment in the ZAZ junction 
is metallic (ACT = 11, class III aGNR). (a) TB-NEGF con¬ 
ductance for the ZAZ junction, (b) DFP conductance repro¬ 
ducing the TB-NEGF result in (a) for the ZAZ junction. The 
central armchair segment controls the Fabry-Perot patterns. 
According to the continuum DFP analysis, the physics un¬ 
derlying such patterns is that of a massless relativistic Dirac- 
Weyl fermionic carrier performing multiple reflections within 
a cavity defined by a double-mass barrier [see diagram in (II)]. 
but with the additional feature that Vi = —Miv% are also 
considered for the leads with zigzag edge terminations (see 
text for details). The mass and Vi parameters used in the 
DFP calculations were L\ — 66ao, mi = 0, Vi = 0, L 2 = lao, 
m 2 V% = 0.50t, V 2 = 0, MiVp = 2 1. Vi = —Miv 2 F . E is the 
energy in units of t. ao = 0.246 nm is the graphene lattice 
constant; t = 2.7 eV is the hopping parameter. 


and mixed armchair/zigzag junctions described in this 
section, because they cannot distinguish between arm¬ 
chair and zigzag edges. 

Discussion 


To motivate the ID Dirac formalism developed and uti¬ 
lized here, we note that customarily armchair or zigzag 
edge terminations in two-dimensional (2D) graphene 
nanoribbons are treated with the massless 2D Dirac 
equation with the use of the corresponding boundary 
conditions.^ In our study, however, narrow GNRs are 
treated with a ID generalized Dirac equation. This ap¬ 
proach is reminiscent of the ID description of GNRs as 
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described in Ref. m and also reviewed in Ref. CD where 
the tight-binding 2D spectrum is projected onto the lon¬ 
gitudinal wavenumber direction, k x . For the armchair 
case (see Eqs. 10 on Ref. ED, the low-energy range of 
the ID spectrum can be well approximated by the Ein¬ 
stein energy relation with a non-zero mass term for the 
semiconducting case, and with a zero mass term for the 
metallic case which exhibits a massless linear dependence 
(photon-like dispersion) of the energy on the momentum. 

As a function of k x . the zigzag GNRs exhibit a partially 
flat E ~ 0 band due to the localized (in the transverse 
direction) edge state. The flat part is centered around 
k x = 7T (alternatively taken as k x = 0 in Ref. 130) , where 
E{k x = 7r) = 0, and it expands towards the graphene 
Dirac points located at 2-7 t/ 3 (alternatively, k x = —7r/3) 
and 47 t/3 (alternatively k x = 7r/3) as the width of the rib¬ 
bon increases (see Fig. 7 in Ref. EJ- After reaching the 
ends of the flat part, the energy band opens two branches 
with non-zero energies; these branches do not have a lin¬ 
ear dependence on k x . Although this can be termed also 
as a “gapless” spectrum (like the metallic armchair case), 
the two differ in their dispersion relation - that is the en¬ 
ergy of the former (zigzag) is approximately independent 
of the momentum for a broad range of the k x momentum 
(see examples for various widths of zigzag GNRs in Fig. 
4 of Ref. [7J , whereas the latter (armchair) shows a linear 
energy-momentum relationship (see Fig. 3 in Ref. E). 

As the width of the GNRs get smaller (narrow zigzag 
nanoribbons), the range of the aforementioned flat part 
decreases, and in the limit of a single row of benzene rings 
(the narrow-most zigzag GNR, referred to as polyacene) 
the dispersion curve transforms® into two (electron-hole) 
parabolas touching at k x = i r. The spectrum of poly¬ 
acene is gapless, but due to the nonrelativistic parabolic 
dispersion, a mass, A4, is associated® 2 with the sec¬ 
ond derivative according to the nonrelativistic relation 
p 2 / 2A4. The case of narrow zigzag GNRs that we study 
here is closer to the nonrelativistic polyacene than to 
the relativistic massless graphene (i.e., the limit of a 
zigzag ribbon of infinite width). Incidentally, we mention 
that the polyacene parabolic-banc! limit is not obtained^ 
through a treatment employing the massless 2D Dirac 
equation with boundary conditions. 

We comment here that physical circumstances where 
a change in some system configuration can result in a 
change in the nature of the dispersion relation (e.g from 
linear to quadratic), are not that rare. Another ex¬ 
ample is found when comparing the gapless, but lin¬ 
ear in momentum, dispersion of the energy of particles 
(massless relativistic) in a monolayer 2D graphene sheet 
with the gapless, but parabolic, dispersion (massive non¬ 
relativistic particles) in bilayer graphene. 34 

Furthermore, our NEGF conductance calculations for 
segmented zigzag GNRs exhibited Fabry-Perot oscilla¬ 
tory patterns, with the spacing between peaks behaving 
as 1 /L 2 (where L is the length of the GNR segment). 
This differs from the (photon-like) FP pattern found 
for the metallic case of armchair GNRs where the spac¬ 


ing between peaks in the oscillatory conductance varies 
as 1/L (characteristic of massless particles with linear 
energy-momentum dispersion, e.g. photons), which were 
studied® by Fabry and Perot. This directly suggests that 
the zigzag GNRs can be described by a nonrelativistic 
limit of the Dirac equation with a sufficiently large mass 
(i.e. by the Schrodinger equation). 

Focusing on the intrinsic properties of the graphene 
lattice, originating from the topology of the honeycomb 
network, and using a NEGF approach, we have stud¬ 
ied here the transport properties of atomically precise 
segmented armchair and zigzag graphene nanoribbons, 
with segments of different widths. Mixed armchair-zigzag 
junctions (with segments of different widths) have also 
been discussed. 

The electronic conductance is found to exhibit Fabry- 
Perot oscillations, or resonant tunneling, associated with 
partial confinement and formation of a quantum box (res¬ 
onant cavity) in the junction. The Fabry-Perot oscil¬ 
lations occur for junctions that are strongly coupled to 
the leads (open system), whereas the resonant-tunneling 
spikes appear for weak lead-junction coupling (closed sys¬ 
tem). In particular, with regard to the FP interference 
patterns, three distinct categories were identified (see the 
Introduction), with only one of them having the char¬ 
acteristics of the optical 21 FP pattern corresponding to 
massless graphene electrons exhibiting equal spacings be¬ 
tween neighboring peaks. 

Perfect quantized-conductance flat steps were found 
only for uniform GNRs. In the absence of extraneous 
factors, like disorder, in our theoretical model, the de¬ 
viations from the perfect quantized-conductance steps 
were unexpected. However, this aforementioned behav¬ 
ior obtained through TB-NEGF calculations is well ac¬ 
counted for by a ID contimuum fermionic Dirac-Fabry- 
Perot interference theory (see Methods). This approach 
employs an effective position-dependent mass term in the 
Dirac Hamiltonian to absorb the finite-width (valence- 
to-conduction) gap in armchair nanoribbon segments, as 
well as the barriers at the interfaces between nanoribbon 
segments forming a junction. For zigzag nanoribbon seg¬ 
ments the mass term in the Dirac equation reflects the 
nonrelativistic Schrodinger-type behavior of the excita¬ 
tions. We emphasize that the mass in zigzag-terminated 
GNR segments is much larger than the mass in semi¬ 
conducting armchair-terminated GNR segments. Fur¬ 
thermore in the zigzag GNR segments (which are al¬ 
ways characterized by a vanishing valence-to-conduction 
energy gap), the mass corresponds simply to the car¬ 
rier mass. In the armchair GNR segments, the carrier 
mass endows (in addition) the segment with a valence- 
to-conduction energy gap, according to Einstein’s rela¬ 
tivistic energy relation [see Eq. @1- 

We observe here that the Dirac Fabry-Perot masses 
that we find to yield agreement with the TB conduc¬ 
tance spectra agree well with those obtained through 
Density Functional Theory (DFT) calculation^ where 
the energy gap A/2 = A ivp for armchair GNRs versus 
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width has been determined. For example, for Af w = 5 
and J\f w = 7, Eq. (1) in Ref. [4] yields A/2 = 0.061 and 
0.281, respectively. These DFT values agree well with the 
DFP values of 0.051 (central segment) and 0.281 (leads) 
given in the caption of Fig. [4j where the H-passivation 
effect according to the DFT was incorporated in the TB 
description. 

The above findings point to a most fundamental under¬ 
lying physics, namely that the topology of disruptions of 
the regular honeycomb lattice (e.g., variable width seg¬ 
ments, corners, edges) generate a scalar-potential field 
(position-dependent mass, identified also as a Higgs-type 
fiel d 28 * 35 (), which when integrated into a generalized Dirac 
equation for the electrons provides a unifying framework 
for the analysis of transport processes through graphene 
constrictions and segmented junctions. 

With growing activities and further improvements in 
the areas of bottom-up fabrication and manipulation 
of atomically precise®^ graphene nanostructures and 
the anticipated measurement of conductance through 
them, the above findings could serve as impetus and 
implements aiding the design and interpretation of 
future experiments. 

Methods 

Dirac-Fabry-Perot model. The energy of a rela¬ 
tivistic fermion (with one-dimensional momentum p x ) is 
given by the Einstein relation E = ^/( p x VF ) 2 + (Mvf) 2 , 
where AI is the rest mass and vf is the Fermi velocity of 
graphene. (In a uniform armchair graphene nanoribbon, the 
mass parameter is related to the particle-hole energy gap, A, 
as At = A/(2vf).) Following the relativisitic quantum-field 
Lagrangian formalism, the mass AI is replaced by a position- 
dependent scalar Higgs field (f)(x) = m(x)vf, to which the 
relativistic fermionic field 4(*) couples through the Yukawa 
Lagrangiarl 35 * Cy = —(/? being a Pauli matrix). For 
4>{x) = 4> o (constant) A ivf = rj> o, and the massive fermion 
Dirac theory is recovered. The Dirac equation is generalized 
as (here we do not consider applied electric or magnetic 
fields) 


[E— V(x)]'i!+ihvFCe— - fUcptx)^ = 0. ( 8 ) 

ox 

In one dimension, the fermion field is a two-component spinor 
4/ = (V , u,V , i) T i u an< A ^ stand, respectively, for the upper 
and lower component and a and /3 can be any two of the 
three Pauli matrices. Note that the Higgs field enters in 
the last term of Eq. | 8 |. V(x) in the first term is the usual 
electrostatic pote ntial, which is inoperative due to the Klein 
phenomenorP®23 and thus is set to zero for the case of the 
armchair nanoribbons (where the excitations are relativistic). 
The generalized Dirac Eq. Q is used in the construction 
of the transfer matrices of the Dirac-Fabry-Perot model de¬ 
scribed below. 

The building block of the DFP model is a 2x2 wave- 
function matrix fl formed by the components of two indepen¬ 
dent spinor solutions (at a point x) of the onedimensional, 
first-order generalized Dirac equation [see Eq. (3) in the main 
paper]. fi play^* the role of the Wronskian matrix W used 
in the second-order nonrelativistic Kronig-Penney model. Fol¬ 
lowing Ref. 1351 and generalizing to N regions, we use the sim¬ 


ple form of fl in the Dirac representation (a = oi, /3 = o 3 ), 
namely 

( e iKx e~ iKx \ 

^k(x') ( _ Ae~^^ x J ’ (^) 

where 


^2 (E - V) - m 2 vf 

“ 2 9 ’ 

n vf 


A = 


TivfK 

E - V + mvf ’ 


( 10 ) 


The transfer matrix for a given region (extending between 
two matching points xi and xi specifying the potential steps 
m( n) ) is the product Mif(a:i,£ 2 ) = fix(x 2 )flf 1 (xi); this lat¬ 
ter matrix depends only on the width X 2 — xi of the region, 
and not separately on xi or * 2 - 

The transfer matrix corresponding to a series of N regions 
can be formecP 5 as the product 


tl,iV+l = W-KiiXiiXi+l), (11) 

i=l,N 


where |a;i+i — Xi\ = Li is the width of the ith region [with 
(m,V,K, A) —> (mi,Vi, Ki, Ai)]. The transfer matrix asso¬ 
ciated with the transmission of a free fermion of mass AI 
(incoming from the right) through the multiple mass barriers 
is the product 

T(E) = n t : 1 (x N+1 )t 1 ,N+i^lk(xi), ( 12 ) 


with k = y/(E - V) 2 - M 2 vj,/{hv F ), | E - V\ > Mvf-, 
for armchair leads V = 0, while for zigzag leads V = 
^Mvf. Naturally, in the case of metallic armchair leads, 
k = E/(Tivf), since AI = 0. 

Then the transmission coefficient T is 


T = 



(13) 


while the reflection coefficient is given by 


R = 



(14) 


At zero temperature, the conductance is given _ by 

G = ( 2e 2 /h)T ; T is the transmission coefficient in Eq. (131. 


TB-NEGF formalism. To describe the properties of 
graphene nanostructures in the tight-binding approximation, 
we use the hamiltonian 

Htb = —t c\cj + h.c ., (15) 

<i,i> 


with <> indicating summation over the nearest-neighbor sites 
i,j. t = 2.7 eV is the hopping parameter of two-dimensional 
graphene. 

To calculate the TB-NEGF transmission coefficients, the 
Hamiltonian ( | 15| ) is employed in conjunction with the well 
known transport formalism which is based on the nonequilib¬ 
rium Green’s functions.®* 

According to the Landauer theory, the linear conductance 
is G(E) = (2 e 2 /h)T(E), where the transmission coefficient 
is calculated as T(E) = Tr[r lQLrQ^]. The Green’s function 
G{E) is given by 

g(E) = (E + i V -Hf%-X l-Xr)- 1 , (16) 
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with H^b being the Hamiltonian of the isolated device (junc¬ 
tion without the leads). The self-energies ^l(r) are given 
by = t l( r)[E + ir\ - tfTB K) ] _lr i(H)> where the hop¬ 


ping matrices tl(r) describe the left (right) device-to-lead 
coupling, and H T b ; is the Hamiltonian of the semi-infinite 
left (right) lead. The broadening matrices are given by 
T L (r) = i[Z L ( R ) - £^ (iJ) ]. 
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